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ABSTRACT 

The Local Group galaxy Phoenix has the properties of a dwarf spheroidal galaxy, but an adjacent 
HI cloud has been recently found to be at the same radial velocity as the stars. The proximity suggests 
that this cloud is associated with the most recent (< 100 Myr) star formation in Phoenix. We have 
obtained relatively high sensitivity and high resolution HI imaging with the VLA with the goal of 
distinguishing between different processes for displacing the gas from the galaxy. Due to the outer 
curvature of the HI cloud, it appears that expulsion from the galaxy by winds from supernovae is more 
likely than ram-pressure stripping. The isolation of the galaxy makes tidal stripping highly unlikely 
Using a star formation history constructed from HST imaging, we construct a simple kinematic model 
which implies that the HI cloud is still gravitationally bound to the galaxy. Gas which is expelled 
from the centers of dwarf galaxies but which remains gravitationally bound may explain the episodic 
star formation observed in several dwarfs. In the specific case of Phoenix, there may be future star 
formation in this currently dSph-like galaxy. 

Subject headings: galaxies: individual (Phoenix dwarf ) — galaxies: dwarf — galaxies: ISM — galaxies: 
kinematics and dynamics — galaxies: evolution 
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1. INTRODUCTION 

The Local Group dwarf galaxies are an intriguing puz- 
zle. Most of them can be divided into two different fam- 
ilies: the dwarf spheroidals (which show no HI gas, no 
current star formation, and no sign of rotational sup- 
port) and the dwarf irregular galaxies (which are HI 
rich, and show obvious signs of recent star formation and 
a significant degree of rotational support). In the last 
decade it has become clear that the dwarf spheroidals 
show a large range in star formation histories — some 
of them are dominated by old (~10 Gyr) stellar popu- 
lations, but many have large components of intermedi- 
ate age (^5 Gyr) and even some young (<1 G yr) stars 
(jMateol 119981 : lGrebelll20?)l iDolnhin et al.il2006D . Stellar 
abundance patterns also show that in many cases the 
process of recycling interstellar gas int o new stars must 
have c ontinued oyer mul t iple Gyr (e.g. , Ikuta k Arimotol 
[2001 : iVenn et al.l l20f)i iFenner et all 12006^ Thus, 5 
billion years ago, many of the galaxies which we call 
dwarf spheroidals must have been significantly more gas- 
rich than they are now and they may have looked very 
much like today's dwarf irregulars. What mechanism or 
mechanisms are responsible for transforming the dwarf 
spheroidals into their current state? 

The fate of the neutral gas in the smallest galaxies 
is relevant to cosmological problems of current interest, 
such as the chemical enrichment history of the early uni- 
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verse. Additionally, it may offer insight into the "dwarf 
galaxy number crisis" - the prediction from standard cos- 
mologies that hundreds of unidentifi ed dark matter ha- 
los should inhabit the Local Group (jMoore et al.l 119981 
iKlypin e~ al. 1999). A long favored theory holds that 
bursts of star formation can remove the ISM from a 
dwarf galaxy be cause of the small gravitational poten- 
tial d epth (e.g., lLarsonlll974t I Vaderl I1989 : iDekel k SiUd 
1986). However, recent hydrodynamical simulations of 
bursts in dwarf spheroidal systems imply that removing 
all of the ISM may be diffi c ult to accomplish in p ractice 
(|Mac Low k Ferraralfl999t iMarcolini et all 12009 ). This 
theory may also require some delicate fine-tuning in order 
to explain a galaxy like the Carina dwarf, which shows 
evid ence of three well isolated episodes of star forma- 
tion (IS aha et al.lll989:lMighel]||1990tlSmecker-Hane et all 
11994 iMighellll 1997T) If the star formation activity is en- 
ergetic enough to clear gas from the galaxy, then cither 
fresh gas was re-accreted between episodes or the clearing 
is only partial. 

Ram pressure stripping presents a second viable the- 
ory, and it has the advantage of naturally explaining 
the morphology-density relationship observed for the Lo- 
cal G r oup dwarfs fe.g.. lEinasto et al.lll974l: iLin fc Faberl 
119831: iKormendvl [19851: Ivan den Berghl U994aUbD. ReT 
centlv. iMaver et all (|2001allbl . I2006D have suggested that 
ram pressure plus "tidal stirring" (the combined effects 
on a satellite galaxy during passage through perigalacti- 
con) can convert a dwarf irregular into a dwarf spheroidal 
by inducing a burst of star formation while simultanc- 
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ously heating the stellar disk and pulling off the ISM 
in a tidal stream. These processes are most important, 
though, on orbits with pericenters less than 50 kpc. 

High resolution and high sensitivity HI observations 
provide a probe to examine which processes are respon- 
sible for removing the gas from the galaxy. For ram pres- 
sure stripping, one expects the gas to be displaced in the 
opposite direction to the galaxy's motion. One might 
also expect that on the side of the cloud facing the op- 
tical galaxy there would be a stronger column density 
gradient giving a bow shoc k appearance si milar to Ho II 
(iBureau fe CarignanH"2002ft and the LMC (|de Boer et all 
1998). For blow-away, naively, one expects gas ejected 
in more than one direction and possibly the presence of 
a bow shock pointed away from the galaxy. For tidal 
stirring, one expects an elongated tidal feature. 

The Phoenix dwarf offers a unique opportunity in the 
study of the mechanisms that can transform a gas-rich 
dlrr into a gas-poor dSph. Its stellar population is pri- 
marily old, but some star for mation activity has contin- 
ued until about 100 Myr ago (lOrtolani fe Grattonlll988l : 
IMartmez-Delgado. Gallart. fc Apariciol 1 1999ft . There is 
also an HI cloud which is associated with the galaxy 
(lYoung fc Ldll997tlsFGermain et ail 1999t [Gallart et all 
12001ft but displaced from the stellar body, and this ar- 
rangement suggests that whatever processes are respon- 
sible for gas removal might be taking place right now in 
Phoenix. The other nearby dwarf galaxies have either 
lost their gas long ago, so that no signs of the removal 
process remain, or they still retain their gas, so that no 
signs of the removal process are visible yet. 

The first map of the —23 km s _1 HI cloud now believed 
to be associated with the Phoenix dwarf was presented 
bv lYoung fc Lol (fl997ft at 133" x 102" resolution. Sub- 
sequent mosaic imaging bv lSt-Germain et al.l (|1999f ) had 
comparable sensitivity and resolution (180" x 180") but 
much better recovery of large scale (> 10') HI structures 
in the region. At the time, the radial velocity of the 
Phoenix dwarf was n ot known, but more recent velocities 
from stellar spec tra ([Gallart et alJl200lt llrwin fc Tolstovi 
120021 lHeldll2002ft showed an excellent match with the HI 
radial velocity, securing the association. Comparisons of 
the properties of the —23 km s _1 cloud with the other 
nearby gas clouds provided evidence that the —23 km s _1 
cloud is indeed physically associated with the Phoenix 
dwarf and that the others are almost certainly foreground 
emission. We now present very deep VLA integrations on 
the —23 km s _1 cloud which give a factor of three better 
sensitivity at improved resolution (71" x 71"). We use 
these high quality images to study the detailed structure 
of the cloud and to seek clues, as described above, to the 
physical mechanisms responsible for the evolution of the 
interstellar medium in the Phoenix dwarf. 

The distance to the Phoenix dwarf has been very 
well constrained by observations of its tip of the 
red giant branch s tars s ince the pioneering work of 
Ivan de Rvdt et all (|1991l ). who first derived a dis- 
tance modulus of 23.1 ± 0.1. Subsequent stud- 
ies have all been consistent with this value as 
IMartmez-Delgado. Gallart. fc Apariciol (|1999ft derived 
23.0 ± O.l. lHeld et alJ ([1999ft derived 23.04 ± 0.07 (and 
a value of 23.21 ± 0.08 base d on the lumino s ity of the 
horizontal branch stars) , and iHoltzman et al.l (|2000ft de- 
rived a value of 23.11 with an uncertainty of "several 



percent." As described in §4, we have re-analyzed imag- 
ing from the HST archive, and find a distance modulus 
of 23.11 ± 0.05 (or, roughly 420 ± 10 kpc). We will use 
this distance throughout the paper. 

2. OBSERVATIONS AND DATA REDUCTION 

The first interferometric HI observations of the 
Phoenix dwarf consisted of 2.0 hours with the National 
Radio Astronomy Observatory's Very Large Array 1 in 
the DnC configuration on 1996 May 18 and 4.0 hours in 
the CnB configuration on 1996 Jan 23. The results of 
those data are published in [Young fc Lol (|1997[ K 

Additional observations in 2003 and 2004 were, in to- 
tal, four times longer than the 1996 data. The Phoenix 
dwarf was observed for a total of 12.0 hours in the DnC 
configuration on 2003 Jan 24, 25, and 27 and for an- 
other 12.0 hours in the CnB configuration on 2004 Jan 
23, 24, and 25. All data were obtained with a total 
bandwidth of 1.56 MHz, giving 256 channels of 6.1 kHz 
(1.3 km s _1 ) each. The velocity range covered is +188.7 
km s _1 to —74.2 km s _1 (heliocentric, in the optical def- 
inition). Integration times were 60 seconds on source in 
the DnC configuration and 20 seconds on source in CnB. 
The nearby point source J0155-408 was observed every 40 
minutes as a phase reference; the bandpass and flux cal- 
ibration data were obtained from the source J0137+331. 

All data calibration and image formation were done 
using standard calibration tasks in the AIPS package. 
Initial imaging revealed which channel ranges were free 
of HI line emission. Continuum emission was subtracted 
directly from the raw visibility data by making first or- 
der fits to the line-free channels. All datasets were then 
combined. For a low resolution map (described in greater 
detail below), the datasets were combined with no change 
to their nominal weights. For a higher resolution map the 
data weights in the 2004 CnB observations were multi- 
plied up by a factor of two. This upweighting partly 
offsets the lower weights resulting from the shorter inte- 
gration times in the CnB array so that after final gridding 
in the visibility plane, both DnC and CnB configurations 
contribute roughly equally to the final image. 

The primary imaging difficulty for VLA imaging of the 
Phoenix dwarf is in achieving a nearly round synthesized 
beam for a source of such low declination. The beam 
is elongated in the north-south direction even for data 
from the DnC configuration, so, in effect, the north-south 
baselines from the CnB configuration are used in combi- 
nation with the east-west baselines from the DnC con- 
figuration to create a synthesized beam which is nearly 
round. A low resolution cube was made using the "nat- 
ural" weighting scheme (data weights do not depend on 
the local density of samples) applied with a Gaussian ta- 
per to downweight long baselines. The resulting beam 
had a fitted FWHM of 130" x 116". The dirty image 
was cleaned down to a residual level of approximately 
1.0 times the rms noise fluctuations, and clean compo- 
nents were restored with a circular beam of FWHM 130". 
A higher resolution cube was made using Briggs's ro bust 
weighting scheme ([Briggs. Schwab, fc Srameklll999f ). as 
implemented in the AIPS task IMAGR with a robust- 
ness parameter 2, plus an additional Gaussian taper to 

1 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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TABLE 1 
Image Parameters 





high res 


low res 


uvtapers (kA) 


3, 50 


1, 3 


robustness 


2 


5 


"original" beam FWHM (") 


71 x 52 


130 x 116 


"original" beam PA (°) 


14 


57 


restoring beam FWHM (") 


71 x 71 


130 x 130 


Linear resolution (pc) 


140 


260 


rms noise (mjyb -1 ) 


1.2 


1.9 


N(HI) sensitivity (cm -2 ) 


2.0X10 18 


9.6xl0 17 



Note. — Column density sensitivities are 3<r in one 
channel. 



downweight baselines extended in the east- west direction. 
The synthesized beam for this image had a fitted FWHM 
of 71" x 52". It was cleaned in the same manner as the 
low resolution cube and restored with a round beam of 
FWHM 71". Offline Hanning smoothing produced chan- 
nels of width 2.6 km s -1 for both cubes. Additional 
imaging details and final sensitivities are given in Table 

m 

Integrated intensity (zeroth moment) maps were made 
by the masking method: the deconvolved image cube 
was smoothed along both spatial and velocity axes, and 
the smoothed cube was clipped at about 1.0a (the noise 
level before smoothing) in absolute value. The clipped 
version of the cube was used as a mask to define a 
three-dimensional volume in which the emission is in- 
tegrated over velocity. Integrating the "moment zero" 
image again over the spatial directions then gives the 
total HI flux. 

The low resolution HI cube contains a total of 2.95 
±0.05 Jykms -1 of HI emission in the vicinity of the 
Phoenix dwarf. The uncertainty in this value is indica- 
tive of the variability that results from slightly different 
sum regions. The higher resolution HI cube contains a 
total of 2.8 ±0.1 Jykms -1 , which is entirely consistent 
with the flux in the low resolution cube. Both values 
are just a bit highe r than the 2.6 Jykms -1 measured by 
lYoung fc Lol Jl997) in their lower sensitivity data cubes. 
Accounting for uncertainty in the absolute flux calibra- 
tion, we count a total flux of 2.95 ± 0. 10 Jykms -1 from 
the P hoenix dwarf. In comparison, ISt-Germain et all 
(1999) quote 4.0 Jykms -1 in this cloud, so the new VLA 
maps have recovered at least 75% of the HI flux. At 420 
kpc, this flux corresponds to 1.2xl0 5 M Q of HI. 

3. THE SPATIAL AND KINEMATIC STRUCTURE IN THE 
PHOENIX HI CLOUD 

As in previous observations of the region, Galactic 
emission is detected in the velocity rang e ±16 to —2 
km s - 1 . The southern clouds no ted by lYoung &: Lol 
(|1997h and ISt-Germain et all £l999) are also evident in 
the velocity range ±44 to ±75 km s -1 . However, since 
these clouds are not believed to be associated with the 
Phoenix dwarf we do not discuss them further. 

Our new images confirm the general impressions gained 
from earlier, lower resolution and lower sensitivity data 
but provide important new details on the internal struc- 
ture of the gas cloud associated with Phoenix. For ex- 
ample, Figure [1] shows the integrated HI intensity of the 
—23 km s -1 cloud at 71" resolution. At column densities 



~10 19 cm -2 the angular exte nt of the cloud is roughl y 
8'x 8', similar to that noted bv I St- Germain et aTl |l999). 
However, the curved shape of the cloud is much more 
obvious here than in the low resolution images. The 
new data also reveal that the "center of curvature" of 
the cloud is not coincident with the optical center of the 
galaxy; rather, it is some 2' to 3' (260 to 390 pc) south 
of the optical center. However, the northern-most and 
southern-most extensions of the HI cloud do lie perfectly 
symmetrically above and below the declination of the 
center of the optical galaxy. Note that the small appar- 
ent optical diameter of Phoenix is a bit deceptive. The 
stellar distribution has be en traced out to a radius of 8.7' 
(|van de Rvdt et ail 1 19911 ). so that the HI observed here 
is within the projected area of the Phoenix galaxy. 

Figure [2] shows individual channel maps at 71" resolu- 
tion; the curved shape of the HI in many of the channels 
and a gradual north-south velocity gradient are clearly 
evident. The peak brightness temperature (1.5 K) is not 
in the center of the cloud but in its southwest quadrant 
near 0l' 1 50™ 44 s , -44° 29' 30". The peak column den- 
sity (5.4x 10 19 cm -2 , or 0.6 M Q pc -2 including He) is also 
here, so that the steepest radial gradients in column den- 
sity are on the side opposite the stellar body. The peak 
brightness temperature at 71" resolution is only modestly 
higher than the brightness temperature at 133" x 102" 
resolution (1.4 K) even though the beam area is smaller 
by a factor of nearly 3, which suggests that the primary 
structures in the HI cloud are now resolved. 

Figure [3] shows the velocity field of the HI emission; it 
also confirms general feat ures noted in the channel maps 
and in the velocity field of lSt-Germain et all (|1999f ). The 
HI cloud shows a velocity gradient primarily in the north- 
south direction from —18 km s -1 on the north side to- 
wards —28 km s -1 in the south. Beyond —23 km s -1 the 
isovelocity contours bend towards the southeast, so that 
the most negative velocities are actually in the south- 
east corner of the cloud (east of 01 51 m and south of 
—44° 29'). The integrated intensity image (Figure [1} 
also shows a protuberance on the west side of the HI 
cloud (01 h 50 m 30 s , -44° 29.5') which departs from the 
velocity trend noted above. Figure [3] shows that the pro- 
tuberance has velocities of —20 km s -1 and higher; it 
is also visible in channel images at —12, —15, and —20 
km s . Even though this feature has quite low col- 
umn density, only lxlO 19 cm -2 at 71" resolution, it 
is clearly detecte d in b oth the ATCA observations of 
I St- Germain et all (|1999f l and in our VLA data, so it is 
undoubtedly real. 

Fits of a Gaussian line profile to the data in velocities 
-2.0 km s -1 to -38.1 km s -1 give the velocity disper- 
sions shown in Figured) Those HI dispersions range from 
8 km s -1 to 14 km s -1 , and, to first order, most of the 
HI is consistent with a value of 10 km s -1 . There is a 
trend for larger dispersions to be found in the northern 
part of the cloud, and, aside from noisy regions at the 
edge of the galaxy, the dispersions are somewhat larger 
on the east side (~ 11 km s -1 ) than on the west (~ 9 
km s -1 ). This trend, with larger dispersions on the side 
closer to the optical galaxy, is consistent with the pattern 
expected from an expanding shell. Looking through the 
outside edge of the shell, the expansion velocity would be 
tangential to the line of sight and would not contribute 
to the line width; closer to the center of the expanding 
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Fig. 1. — Integrated HI intensity at 71" resolution overlaid on the optical image of Phoenix from the DSS. Contours are (—2, —1, —0.5, 
0.5, 1, 2, 3, 4, and 5) xlO 19 cm -2 . The small hashed circle in the top left corner indicates the resolution of the HI image. The crossed 
circle marks the bright part of the stellar distribution and is intende d to guide the eye for a comparison with Figure [2] Note that the 
stellar distribution of Phoenix has been traced out to a radius of 8.7' (van dc Rydt ct al. 1991), as shown by the large circle, so that the 
HI observed here is within the projected area of the Phoenix galaxy. 



shell (on the north and east sides of the HI cloud, in this 
case) the expansion velocity would have some observable 
component and would broaden the line. However, if the 
HI cloud is an expanding shell its expansion velocity is 
only a few km s~ x and it is not possible to distinguish the 
expansion signature in position- velocity slices dominated 
by the overall N-S velocity gradient. 

Neither our presen t data nor the results of 
ISt-Germain et al.l |l999) give any hints for the presence 
of additional HI emission offset in any other direction 
from the optical galaxy, down to column density limits 
of approximately 10 18 cm -2 at 130" resolution (40 times 
smaller than the observed column densities). The signif- 
icance of this non-detection lies mainly in the question of 
the importance of ram pressure as a gas removal mech- 
anism. HI gas offset in multiple directions would almost 
certainly rule out ram pressure stripping as a credible 



explanation; however, we find no such gas. 

In considering the possible effects of ram pressure 
on the gas in Phoenix, though, the morphology of the 
HI cloud is particularly interesting when compared to 
the morphology of HI in t he dwarf galaxy Holmberg II 
(|Bureau fc Cari gnan 2002h. The low resolution map of 
Ho II clearly shows smooth, convex, and compressed HI 
contours on the southeast side of the galaxy, with more ir- 
regular and fluffy low column density material extending 
to larger radii on the north and west sides. The cometary 
appearance of the HI in Ho II is clear evidence for the 
effects of ram pressure as the galaxy moves through the 
M81 intra-group medium. If ram pressure were responsi- 
ble for displacing the HI from the center of Phoenix, we 
would expect to find a steeper HI column density gradi- 
ent on the side towards the stellar body. We might also 
expect the cloud to have a convex side towards the stellar 
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Fig. 2. — Channel maps at 71" resolution. Contours are —3, 3, 6, and 9 times the rms noise (1.2 mjyb" 1 ). The resolution is indicated 
in the top left corner of the first panel, and the crossed circle shows the position and approximate size of the stellar body visible in the 
DSS image of Figure [T] 



body. Neither of those is true, and thus the comparison 
with Ho II suggests that ram pressure is not responsible 
for the displacement of the HI in Phoenix. 

Quantitatively, it is also unlikely that the density of 
the halo gas is high enough to effect any ram pres- 
sure stripping at Phoenix. In the usual manner we can 
estimate the neces sary density from the relation (e.g. 
iGunn fc Gottj fl972T) P ram ~ p ha i v 2 gal ~ GM gal Z gas /R 2 , 
where v ga i is the velocity of the galaxy through the halo 
gas, Mg a i is the mass of the galaxy interior to some ra- 
dius i?, and T, gas is the column density of the atomic 
gas there. Making the most conservative estimates for 
the HI properties of Phoenix we consider the outer edge 



of the HI cloud, which has surface density 0.1 M Q pc 2 
(including helium) at 590 pc from the gala xy center. We 
take Phoenix's mass to be 3.3xl0 7 M (|Mateol fl998f) 
and its velocity about the Milky Way approximately 220 
km s . Stripping then requires the density of the halo 
gas to be phaio ^ 3xl0~ 5 cm" 3 . Stripping the gas of 
the highest observed column densities from the center of 
Phoen ix requi r es at least an order of magnitude larger 
Phaio- iMuralil (|2000f ) has estimated that the density of 
the Milky Way halo must be < 10" 5 cm" 3 at the lo- 
cation of the Magellanic Stream, however, and Phoenix 
is some eight times further from the Milky Way than 
that. Of course ram pressure stripping is easier if the 
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Fig. 3. — Velocity field, derived from fitting Gaussian line profiles 
to the image at 130" resolution. Contours range from —30 to —14 
km s _1 in steps of 2 km s _1 . 
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Fig. 4. — Velocity dispersion at 130" resolution. Contours range 
from 9 km s — 1 to 13 km s — 1 in steps of 1 km s — 1 . Most of the HI 
can be characterized as having a velocity dispersion between ~9 
and 11 km s — . A few pixels in which the formal uncertainty in 
the dispersion is greater than 2.5 km s — 1 are not shown. 



mass of Phoenix has been ov erestimated, but by analogy 
with iMarcolini et al.l ([2006) , if the mass of Phoenix is 
very much smaller than we have assumed here it could 
not have retained its ISM long enough to support an ex- 
tended star formation history. 

Tidal stirring is another process which can be respon- 
sible for the removal of gas from a galaxy. However, 
since this process is most rel evant at distances less than 
50 kpc from the Milky Way (jMaver et al.ll2006l ). the rel- 
ative isolation of Phoenix makes this scenario unlikely. 
The radial velocity of —23 km s _1 corresponds to a galac- 
tocentric radial vel ocity of —112 km s" 1 using the for- 
mulae presented in IdeVaucouleurs et al.l (|1991f ). At 112 
km s _1 it has been at least 3 Gyr since Phoenix was close 
enough to the Milky Way for significant tidal interaction. 
If Phoenix is bound to the Milky Way, the semi-major 
axis of its orbit must be at least 210 kpc, giving an or- 
bital period rfTJ) Gyr for a Milky Way mass of 1.4x 10 12 
M Q (jMateolll99ct ). Clearly, interactions with the Milky 
Way may have influenced Phoenix several Gyr ago but 
not more recently than that. 

4. THE RECENT STAR FORMATION HISTORY OF THE 
PHOENIX DWARF 

In an early groun d -based optical study of Phoenix, 
lOrtolani & Grattonl (jl988) noted the' presence of 
a small (~ 10 4 M Q ) population of young (~ 
10 8 yr) blue stars. 2 From deep e r ph otometry, 
iMartinez-Delgado Gallart. fc AparlcToT (|1999h (MGA) 
studied the distribution of these blue stars and found 
evidence of an east-to-west progression across the face of 
Phoenix. MGA estimated the age of these stars to be 
100 Myr, and, fro m this, the relative position of the HI 
cloud detected by lYoung &: Lol (|1997ft . and the assump- 
tion that the HI cloud was associated with the recent 
star formation, estimated an outward velocity for the HI 
cloud of ~ 7 km s . At the time, the optical radial 
v elocity of Pho e nix w as not known. 

iGallart et al.l (|2001l ) obtained velocities for 31 Phoenix 
red giants and determined an optical radial velocity of 
—52 km s . The ^30 km s _1 difference between the op- 
tical radial velocity and the HI radial velocity required a 
re-calculation of the estimated HI cloud velocity. The ra- 
dial velocity difference implied a trajectory nearly along 
the line of sight, because placing the HI cloud back on 
top of the stars 100 Myr ago means the tangential ve- 
locity of the cloud is much less than 30 km s _1 . The 
resulting galaxy-HI cloud separation is 3.5 kpc. At this 
velocity and separation, the HI cloud is not bound to 
Phoenix (the mass estimate falls short by roughly one 
order of m agnitude). 

However , llrwin fc Tolstovl (12002ft obtained a new radial 
velocity of —13 ± 9 km s _1 , which is in exc ellent agree- 
ment with the HI velocity of -23 km s _1 . iHeldl (2002) 
also reported agreement between the optical and HI ra- 
dial velocities within 2-3 km s . Thus, the original 
calculation by MGA wou l d app ear to be more relevant 
than that of lGallart et~aT1 (2001) concerning the gravita- 
tional status of the HI cloud. 

We have revisited the recent star formation history 
of the Phoenix dwarf in the following way. The 

2 The abstract of the paper states 10 7 yr, but the text assigns 
an age of 10 s yr. 
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Fig. 5. — The V versus VT color magnitude diagram from the 
HST WFPC2 observations of the central field of Phoenix (see text). 
The photometry comes from the Local Group Stellar Populations 
Archive (Holtzman et al. 2006). The right axis shows the absolute 
magnitude scale corresponding to the derived distance modulus of 
23.11 ± 0.05. Note the truncation in the main sequence at an 
apparent magnitude of ~ 21.5. This provides an accurate determi- 
nation of the age of the last star formation in Phoenix. 

Hubble Space Telescope WFPC2 has obtained rela- 
tively deep opti cal imaging of Phoenix in two differ- 
ent programs. iHoltzman et alj (|2000f ) report obser- 
vations of a central field in Phoenix. They derive 
a star formation history for this field, and note that 
star formation proceed ed in Phoenix up un t il 100 Myr 
ago, in agreement with [Ortolani &: Grattor] dl988f) and 
iMartmez-Delgado. Gallart. &: Apariciol (|1999l ) . Deeper 
HST WFPC2 observations of a central field which is 
almost coincident with that of Holtzman et al. were 
obtained by Aparicio (HST-GO-8706) in 2001. We 
have retrieved the photometry from these two Hub- 
ble Space Telescope observations of Phoenix from the 
HST/WFPC2 Local G roup Stellar Photometry Archive 
ijHoTt zman et"ai1l2006l ). For reference, the V,V~I color 
magnitude diagram for the deeper integration is shown 
in Figure [5] 

Using the programs of lDolphir] (|2002f ). we have derived 
the star formation histories for both Phoenix central field 
HST observations. In Figure [5] we show the total star 
formation history for the central field derived from the 
deeper integration. This star formation history agrees 
well with that derived from the slightly shallower inte- 
gr ation and the analysis o f those observations presented 
bvlHoTt zman et a l. (2000). There is evidence for a large 
amount of star formation at the earliest ages, a relative 
decrease at intermediate ages (~ 10 - 5 Gyr ago), and 
then increased star formation activity in the last 5 Gyr 
which decreases as the present time is approached. 

In Figure [7| we show the star formation history for the 
last Gyr derived from the deeper integration at the higher 
time resolution afforded by the observations. The age of 
the last star formation event appears to be about 50 Myr 
ago with a duration of roughly 20 Myr. A similar analysis 
of the shallower central field observation gives a similar 
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Fig. 6. — The global star formation history of the Phoenix dwarf 
galaxy derived from HST observations of a central field (see text). 
This star formation history is in generally good agreement with 
the ce ntral field star formation history derived by Holtzman et al. 

result. The age of 50 Myr for the most recent episode 
of star formation is slightly younger than the previous 
reported values of ~100 Myr ago. However, this is the 
first detailed analysis of this event from HST photometry, 
and Figure [5] shows that there is a clear truncation in the 
main sequence stars at an apparent V magnitude of ~ 
21.5, so the age estimate should be secure. Thus, it is 
likely that 50 Myr is a better estimate than 100 Myr for 
the most recent star formation event. Note that this most 
recent star formation episode is comparable in strength 
to the earlier episode at ~200 Myr ago. (Because linear 
time is plotted on the horizontal axis of Figure [TJ the 
total star formation for any time is simply the area under 
the curve for the appropriate time bin.) Also note that 
the "gaps" in star formation at ~300 Myr and ~100 Myr 
ago are of similar duration (~50 Myr). 

Using the age of 50 Myr for the most recent star for- 
mation event, we can calculate simple estimates of the 
HI cloud velocity and determine its gravitational status. 
The angular separation between the peak column den- 
sity of the HI cloud and the optical center of the galaxy 
is 289". At a distance of 420 kpc, this corresponds to 
a projected distance of 590 pc. Under the assumptions 
that the HI cloud is 590 pc from the galaxy and was 
expelled by the most recent episode of enhanced star for- 
mation, we derive a velocity of ^12 km s^ 1 . An upper 
bound on the velocity and separation of the HI cloud can 
be obtained by including the observed radial velocity dif- 
ference between HI and the stellar body (10 km s _1 ); the 
recession speed of the HI cloud could be as large as 16 
km s _1 and its true separation as large as 780 pc. 

Giv en Phoenix's total mass of 3.3xl0 7 M ([Mated 
QUI, the escape velocity for the HI cloud at the dis- 
tance of 590 pc is 21 km s _1 . At 780 pc from Phoenix, 
the escape velocity is 19 km s _1 . The HI cloud's veloc- 
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Fig. 7. — The recent star formation history of the Phoenix dwarf 
galaxy derived from HST observations of a central field (see text). 
The most recent episode of star formation has an age of roughly 
50 Myr with no evidence of star formation since then. 



ity is 3-9 km s _1 lower than the escape velocity, so we 
find that the cloud is still gravitationally bound. It is 
likely that the HI will fall back and may serve as the 
raw material for future episodes of star formation. How- 
ever, we also note that there is an observed velocity dif- 
ference of 12 km s _1 from one side of the cloud to the 
other, meaning that some portion of the cloud may be 
u nbound. Thes e conc lusions are contradictory to those 
of lGallart et al.1 (|2001l ). but this is purely a result of the 
previous authors' use of a relatively large radial velocity 
offset between the HI and the stars. 

For the observed HI mass, multiplied up by a factor 
of 1.3 to account for helium, and the recession velocities 
discussed above, we find that the total mechanical energy 
in the HI cloud is 2.3xl0 50 - 3.9xl0 50 ergs. Integrating 
over the star formation history, we calculate a total mass 
of ^5xl0 3 M© of new stars formed in the recent star 
formation event. For a standard IMF, one expects more 
than 10 stars with masses greater than 10 M formed in 
that event, so, even at low efficiency, there is sufficient 
energy available from the supernovae associated with the 
recent star formation to propel the HI away from the 
center of the galaxy. 

Finally, there remains the question of why neutral gas 
was only ejected to one side of Phoenix. Naively, one 
expects gas blown away by a burst of star formation to 
be ejected in all directions unless the energy source was 
initially off-center with respect to the gas. Could one 
off-center supernova be responsible for the asymmetric 
gas displacement in Phoenix? Assuming a canonical to- 
tal energy for a supernova of 10 51 erg implies an energy 
transfer efficiency in the range of 20%-40%. Note, how- 
ever, that this efficiency does not take geometric effects 
into account — if the blast energy were injected sym- 
metrically, much less than half of the 10 51 erg would be 



available to accelerate a cloud off to one side of the super- 
nova. Thus, our estimates are inconsistent with energy 
input from a single typical supernova. The idea also ne- 
glects the effects of the other ~ 10 supernovae that likely 
occurred. However, the combined effects of those super- 
novae could be responsible for the single-sided expulsion. 
Based on a comparison between the distributions of the 
brightest main sequence stars and the red su pergiants, 
iMartmez-Delgado. Gallart. fc Aparic io (1999) found ev- 
idence for star formation propagating across Phoenix 
from east to west. The star formation wave may act 
as a plow to sweep up the gas and push it to one side of 
the galaxy. 

We have investigated the HST observations to see 
if the detailed star formation history supports this 
interpretation of single-sided expulsion. The ac- 
curacy of the photometry of the main sequence 
stars in the HST photometry is substantially bet- 
ter thjiii_yiB : t_of_yie^ used 
bv IMartmez-Delgado. Gallart. fc Aparicio! (|1999r i. so it 
is possible for us to examine the spatial distribu- 
tion of main sequence stars as a function of lu- 
minosity (or upper age limit). In agreem ent with 
IMartmez-Delgado. Gallart. fc Aparicio! (|1999f ). we find 
the majority of the brightest main sequence stars fall 
in a band stretching across the middle of Phoenix. The 
brightest of these (stars with indicative ages of 25 - 75 
Myr) are clustered in the association on the west side of 
the galaxy. Stars with older ages are more centrally con- 
centrated. We also derived separate star formation his- 
tories for the eastern side of the galaxy and the western 
side of the galaxy. The western side shows evidence for 
strong episodes of star formation at 40 and 180 Myr ago 
preceded by star formation ~ 400 Myr ago. The eastern 
side shows evidence of relatively constant star formation 
at lower levels from 50 to 250 Myr ago and also star for- 
mation at ~ 450 Myr ago. In sum, the data are consistent 
with a general trend for the star formation propagating 
from east to west over the last 200 Myr, but the pat- 
tern is not a clear-cut, monotonic one (i.e., both sides of 
the galaxy participated in the star formation episode 200 
Myr ago). These data are not inconsistent with the HI 
location on the west side of the galaxy. 

5. DISCUSSION 

The fact that the HI cloud appears to be bound to 
Phoenix is intriguing. From the asymmetrical shape of 
the HI cloud and the relative isolation of Phoenix, it 
appears that ram pressure stripping is not a likely cause 
of its displacement. Thus, supernova winds appear the 
most likely explanation for the HI cloud to be driven 
away from the center of the galaxy. However, if the HI 
cloud is bound, then ~ 10 5 M© of HI can fall back into 
the galaxy, perhaps triggering a future episode of star 
formation. 

Indeed, stellar population work has already suggested 
that this kind of a scenario (temporary HI expulsion, 
fall-back, and rejuvenated star forma tion) may haye hap - 
pened in the Carina dSph. Early on. iMould et all (|1982l ) 
noted that the presence of blue stars in the Carina 
dwarf indicated that Carina and possibly other dSphs 
were not single age systems. Further evidence for a 
complicated star f o rmation histor y for Carina gathered 
(jSaha et all 1 19861 iMighelll 119901 : ISmecker-Hane et all 
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11994 iMighelll I1997D until the present day picture of 
Carina consists of three episodes of star formation 
at roughly old, intermediate (7 Gy r ) and young (3 
Gyr) epochs dHurlev-Keller et al.|[l99alRizzi et aDl2003t 
iMonelli et al.ll2003h . Two major questions which have al- 
ways plagued studies of the star formation history of the 
Carina dSph are: (1) what caused the large gaps between 
the episodes of star formation, and (2) what removed the 
neutral gas from Carina such that there is no longer any 
star formation at all? 

On the second question, the final clearing of gas from 
Carina and indeed the other dwarf spheroidals, Phoenix 
does not provide a direct answer because its HI is ap- 
p arently still boun d. The hydrodynamical simulations 

of iMarcolini et ail (|2006f) and others also show that it is 
difficult for star formation activity, by itself, to blow out 
all of the ISM unless the dark matter halos are signifi- 
cantly smaller than expected and/or the hot gas's radia- 
tive energy losses much smaller. The simulations do tend 
to lift a significant component of the galaxy's cold gas up 
away from the stellar body in the aftermath (30 Myr - 50 
Myr) of a burst, as we observe in Phoenix. But if star for- 
mation cannot completely unbind the ISM from a dwarf 
spheroidal, it is commonly assumed that some external 
agent such as tidal forces or ram pressure stripping must 
unbind the gas. This idea is certainly consistent with 
the fact that the most gas-poor dwarf spheroidals such 
as Draco and Ursa Minor tend to be closer to the Milky 
Way or to M31 than Phoenix is. A recent proper mo- 
tion estimate for the Ursa Minor dwa rf spheroidal implie s 
that its perigalacticon is < 50 kpc (|Piatek et al.ll2005T ). 
carrying it well into the region whe re interactions wit h 
the Milky W a y can clean it of its gas (|Maver et al.ll2006f ). 
iPiatek et all (pOSh have made the interesting comment, 
however, that perigalacticon distance cannot be the only 
factor influencing the evolution of a galaxy as the Ca- 
rina and Ursa Minor dwarfs have similar orbits but quite 
different star formation histories. 

With regard to the first question posed above, the ori- 
gin of the gaps in the star formation activity of Carina, 
the Phoenix dwarf may be providing us with some an- 
swers. Here we posit that a relatively modest amount of 
star formation has lifted the neutral gas up out of the cen- 
ter of the potential well. However, most of the gas is not 
moving in excess of the escape velocity and is destined 
to return to the center. Again, from the separations and 
velocities discussed in S|4]we can calculate the total en- 
ergy (gravitational plus mechanical) of the HI cloud, the 
semi-major axis of its orbit and its orbital period. The 
estimated orbital periods are in the range 140 Myr to 660 
Myr, so it seems reasonable that this mechanism could 
produce gaps in the star formation history on timescales 
of a few hundred Myr (Figure [7]) . On the 660 Myr orbit 
the semi-major axis is 1.2 kpc; thus, the HI cloud could 
reach angular separations several times larger than the 
current value. 



Given the present distribution of the HI gas, it seems 
unlikely that it will all fall back in at once, so there may 
be some additional time delay as the gas collects before 
initiating the next episode of star formation. Since only 
~ 10 3 Mq of stars were formed in the episode that dis- 
located the gas, we may expect a similar amount of star 
formation when the ~ 10 5 M Q of HI gas returns to the 
galaxy's center. Thus, the Phoenix dwarf may be aptly 
named. 

6. SUMMARY 

We have presented new VLA HI observations of the Lo- 
cal Group galaxy Phoenix. The galaxy has the properties 
of a dwarf spheroidal galaxy, with the exception that an 
HI cloud is located at a similar velocity and within the 
optical radius of the stellar distribution. The close as- 
sociation of the HI cloud with the galaxy suggests that 
this cloud may be the remnants of gas responsible for the 
most recent star formation in Phoenix which took place 
< 100 Myr ago. 

Our new high sensitivity and high resolution VLA HI 
imaging reveals that the HI distribution has a crescent- 
like shape with the center of curvature in the general 
direction of the galaxy. This indicates that ram pressure 
stripping is unlikely to be responsible for the displace- 
ment of the gas from the galaxy. Additionally, the isola- 
tion of the galaxy makes tidal stripping highly unlikely. 

This leaves expulsion from the galaxy center by winds 
from supernovae as the most likely explanation. Using a 
star formation history constructed from HST imaging, we 
construct a simple kinematic model which implies that 
the HI cloud is still gravitationally bound to the galaxy. 
Much of the atomic gas can be expected to collect again 
in the center of the galaxy on timescales of a few hundred 
Myr. Thus, there may be future star formation in this 
dSph-like galaxy. 
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